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ABSTRACT 


Based  upon  the  classical  vortex  system  representation, 
the  induced  velocities  of  a  finite-bladed  propeller  with 
arbitrary  circulation  distribution  are  derived  and  Fourier 
analyzed.  The  zeroth  harmonic,  or  steady  component,  of  the 
induced  velocities  is  considered  in  detail  and  formulas  involv¬ 
ing  only  an  integration  over  the  blade  radius  are  found.  The 
mathematical  equivalence  between  these  results  and  the  conven¬ 
tional  actuator  disk  representation  of  the  propeller  is  demon¬ 
strated.  Sample  calculations  of  both  the  axial  and  radial 
velocity  components  for  an  arbitrary  representative  circulation 
distribution  closely  approximating  the  Goldstein  optimum  are 
presented.  For  the  special  case  of  a  uniformly  loaded  pro¬ 
peller,  expressions  for  the  induced  velocity  components  are 
given  in  simple  closed  form  and  velocity  profiles  are  compared 
with  the  representative  results. 


2 


INTRODUCTION 

The  problem  of  developing  suitably  simple,  yet  accurate, 
expressions  for  the  induced  velocity  components  of  a  free 
propeller  has  long  occupied  the  attention  of  both  the  aerody- 
namicist  and  the  naval  architect.  Detailed  knowledge  of  the 
velocity  field  is  extremely  important  for  many  applications 
ouch  as  the  design  of  ducted  propeller®  for  VTOL  aircraft,  the 
evaluation  of  the  interaction  between  a  propeller  and  a  wing 
immersed  in  its  slipstream,  the  prediction  of  the  loading  on 
nearby  appendages,  and  the  determination  of  the  free-space 
sound  field. 

Early  approaches  to  the  development  of  airscrew  theory 
followed  two  lines:  one,  the  momentum  theory  of  R.  E.  Froude  [1] 
and  w.  J.  M.  Rankine  [2]  and  the  related  actuator  disk  con¬ 
cept:  and  the  other,  the  blade-element  theory  of  W.  Froude  [3] 
and  S.  Drzewiecki  [4].  While  able  to  predict  propeller 
performance  with  some  reasonable  accuracy,  the  utility  of  each 
theory  is  limited  by  basic  deficiencies.  On  the  one  hand, 
momentum  theory  yields  neither  a  relationship  between  the 
propeller  geometry  and  the  thrust  and  torque  nor  any  detailed 
knowledge  of  the  induced  velocity  field.  On  the  other  hand, 
blade-element  theory  suffers  from  the  uncertainties  involved 
in  assuming  the  aerodynamic  characteristics  of  the  blade 
sections.  It  also  gives  no  information  regarding  the  induced 
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velocity  field. 

Prediction  of  the  induced  velocity  field  and  in  partic¬ 
ular  the  local  inflow  at  the  propeller  plane  was  finally 
achieved  with  the  now-classical  vortex  theory  of  the  propell¬ 
er.  Analyzed  in  detail  by  A.  Betz  and  L.  Prandtl  [5]  and 
S.  Goldstein  [6]  and  refined  by  others  [7]  &  r8],  it  has 
enabled,  in  conjunction  with  blade-element  theory,  the 
accurate  determination  of  propeller  performance  for  the 
forward  flight  regime.  Still,  most  calculations  of  induced 
velocities  using  vortex  theory  are  either  so  complex  as  to 
require  extensive  time  on  high-speed  electronic  computers, 
can  not  be  readily  carried  out  for  arbitrary  distributions  of 
blade  circulation,  or  arc  based  on  the  oversimplifying  assump¬ 
tion  that  the  number  of  blades  is  very  large  and  the  flow 
periodicity  may  bs  neglected. 

The  purpose  of  this  paper  is  to  remove  some  of  the 
undesirable  features  of  these  theories,  the  results  being 
essentially  an  extension  of  part  of  an  earlier  analysis  for 
the  ducted  propeller  by  D.  E.  Ordway,  M.  M.  Sluyter  and 
B.  U.  0.  Sonnerup  [9].  In  outline,  we  consider  a  lightly- 
loaded  propeller  of  arbitrary  blade  number  and  circulation 
distribution  operating  at  zero  incidence  in  a  uniform  free 
stream.  The  propeller  is  represented  by  the  conventional 
vortex  system  and  the  induced  velocities  at  any  field  point 
are  determined  from  the  Biot-Savart  law.  These  velocities 
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are  then  Fourier  analyzed.  The  zeroth  harmonic,  or  steady 
component,  is  considered  in  detail  and  identified  with  the 
actuator  dish,  or  infinite-blade-number  solution.  Simplified 
formulas  involving  only  an  integration  over  the  blade  radius 
are  derived  and  sample  calculations  for  an  arbitrary  represen¬ 
tative  circulation  distribution  are  carried  out.  For  the 
special  case  of  uniform  propeller  loading,  the  integrations 
are  performed  analytically  yielding  simple  clcsed-form  express¬ 
ions  for  the  velocity  field.  Several  profiles  are  calculated 
and  compared  with  those  for  the  representative  case. 


BASIC  FORMULATION 

consider  a  propeller  operating  at  zero  incidence  in  a 
uniform,  inviscid,  incompressible  stream  of  speed  U  .  We 
assume  the  propeller  to  be  lightly  loaded  with  N  equally 
spaced  blades  of  radius  R  which  are  rotating  about  their 
axis  at  a  constant  angular  velocity  fi  .  In  addition,  we 
assume  that  both  the  blade  thickness-to-chord  and  chord-to- 
radius  ratios  are  negligibly  small  and,  for  simplicity,  dis¬ 
regard  the  hub. 

A  cylindrical  propeller-fixed  coordinate  system  (x,r,e) 
is  chosen  such  that  the  propeller  dish  is  normal  to  the 
x-axis  and  is  located  at  x  =  0  ,  see  Fig.  1.  In  accordance 
with  our  assumptions,  v*e  can  new  represent  the  propeller  by 


g^asaasKM 


>T,f»  <  ? 


<■" 

o 

the  classical  model  of  a  bound  radial  vortex  line  of  strength 
r(rv)  for  each  blade  accospanied  by  a  helical  sheet  of  vor¬ 
tices  of  strength  -  dr/dry  trailing  from  each  line,  where 
rv  is  the  radial  distance  to  any  element  of  the  vortex  line. 
Also  consistent  with  our  assumptions,  the  helical  path  of  the 
trailing  vortex  system  is  determined  solely  by  the  incoming 
free  stream  with  translation  U  and  rotation  ft  . 

In  order  to  express  the  induced  velocity  field  it  is 
convenient  to  introduce  three  elementary  vector  velocity 
fields  [9].  These  fields  or  influence  functions  ,  1^. 
and  are  due  to  vortices  of  unit  strength  and  unit  length 
which  lie  in  the  axial,  radial  and  circumferential  directions 
respectively.  From  the  Biot-savart  law  we  have, 

lx,  ...  -  i  X  D  /  M£|3  (1) 

where  _i  is  the  unit  vector  in  the  poritive  direction  of  our 
desired  element  and  D  is  the  vector  from  the  element  at 
(xv,rv,0v)  to  the  field  point  (x,r,e)  .  With  i^  ,  ir  and 
ig  corresponding  to  the  unit  vectors  in  the  x ,  r  and  6 
directions  at  the  field  point,  the  influence  functions  reduce 
to 


-rvsin(e-9v)j^.  +  [r-rvcos(e-ev) ]ig 
4ttD3 
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h 


r  810(0-6^1*  -  (x-x^BlnfS-Gj^.  -  (x-x^coafG-e^ 


With  these  results,  the  total  velocity  field  for  both  the 
bound  and  trailing  vortices  may  be  simply  expressed.  For  Op  or 
the  velocity  induced  by  the  N  bound  blade  vortex  lines,  we 
integrate  over  each  blade  and  sum  the  result  over  all  blades,  or 

N  R 

3 r  =  Z  /  r<rv>  IrtO'V2*1/1*)  ^v  (4) 

J»1  0 


Here,  the  argument  of  is  the  location  of  the  vortex  element, 

the  field  point  remaining  arbitrary  for  the  present,  and  the 
index  i  denotes  the  ordinal  number  of  the  blade. 

For  or  the  velocity  induced  by  the  trailing  vortices, 

ve  firBt  combine  an  axial  and  a  tangential  element  appropriately 
to  form  an  arbitrary  helical  element  lying  along  the  free  stream. 
Then  we  integrate  over  the  vortex  trailing  from  each  radial 
location,  carrying  out  a  second  integration  over  the  radius  and 
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sunsning  over  the  propeller  blades,  we  arrive  at 

NR  » 

%'  "  I  /  -  r'(r„)  /  <°  4  ♦  <*„  4)  4t  drv  (5) 

i=l  0  0 

where  we  have  set  r'  =  dr/drv  and  the  arguments  of  1^  and 

h  are 

=  Ij^UT.r^i/N  +  oT) 

=  le(UT,rv,27ri/N  +  fir)  (6) 

The  dummy  variable  t  represents  the  time  for  a  trailing 
element  convected  by  the  free  stream  to  travel  from  the  blade 
to  its  position  downstream  and  the  corresponding  integral  over 
t  is  equal  to  the  velocity  induced  by  a  semi-infinite  helical 
vortex  of  unit  strength,  radius  ry  and  pitch  U/'ftrv  . 

Superposition  of  the  velocity  fields  given  by  Eqs.  4  anu  5 
determine  then,  within  the  limits  of  our  formulation,  the  total 
induced  velocity  at  any  field  point.  Within  the  propeller 
slipstream  these  integrals  are  to  be  interpreted  in  the  Cauchy 
Principal  value  sense. 

FOURIER  ANALYSIS  OF  VELOCITY  FIELD 


In  the  coordinate  system  translating  but  not  rotating  with 
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the  propeller,  the  flow  is  periodic  with  time.  That  is,  the 
velocity  at  any  field  point  will  consist  of  a  steady  or  time- 
dependent  component  plus  a  superimposed  fluctuating  component 
which  has  a  period  of  2ff/NH  •  Accordingly,  we  can  express 
and  *n  the  complex  Fourier  series 

2p  °  &  )  Sp  (x»r)  e<mN9 

R\ 


V 


00 


Sr. 


(x,r)  e 


imN0 


m=-« 


m 


(7) 


in  our  propeller-fixed  coordinates  in  which  the  angular  varia¬ 
tion  is  equivalent  to  the  time  dependency.  The  complex  Fourier 

velocity  coefficients  c-  and  c^,  are 

~ 1  m 


-7T 


-rr 


(8) 


from  orthogonality,  U  being  introduced  for  dimensional  pur¬ 
poses  . 

At  this  point  it  is  convenient  to  decompose  , 

c_  and  c_^  into  their  axial,  radial  and  tangential  compo- 

~i  m  ~  m 

nents  respectively,  or 


and  the  components  of  the  complex  Fourier  velocity  coeffi¬ 
cients  associated  with  u_  ...  will  be  denoted  by  cJJ  .... 

m 

To  find,  say,  the  axial  component  cn  ,  we  replace  l_ 

rm 

by  I^«i^  in  Eq.  4  using  the  second  of  Eqs.  2  .  From  the 
first  of  Eqs.  9  the  resulting  expression  represents  u^,  .  If, 
in  turn,  we  replace  by  this  u^,  in  the  first  of  Eqs.  8 

and  reverse  the  order  of  the  rv  and  0  integrations,  we  can 
carry  out  the  6  integration  in  terms  of  Legendre  functions 
[10].  Noting  that  each  blade  contributes  equally  to  any 
harmonic  for  N  identical  blades,  we  obtain 


R 

-iN  f 
87r2U  «/r 


'  3/2  [®mN+l/2  W  ^  "  QmN-3/2^1^]  drv 
v 

(10) 


where  0mN+1^2  and  QmN_3/2  are  the  Legendre  functions  of 
the  second  kind  and  half-integer  order  with  argument  , 


=  1  + 


*2+(r-rv)2 

2rr 


(11) 


and  the  prime  ( ' )  denotes  differentiation  with  respect  to 
this  indicated  argument. 
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The  remaining  coefficients  c£  ...  are  obtained  in 

1  m 

similar  fashion  and  are  given  by 

C^m  “  S^ur3^  f  ~37?  [W(<Dl)  -  Q«/2K)]  drv 


0  *v 


f  r  3/2  [QmN+l/2^1^  +  QmN-3/2^<ul)]  drx 

•*  mm 


m  bffTJr- 


C  *v 


B  r.' 


_u  _  -NfJ 

Cr; 


i»  p  /  ^  j  00 

I  \rzr-  f  frQmN+l/2^m2^  +  rQmN-3/2^m2^ 

0  vrv  0 

"  2rv^l/2(^]e"tmSir  dT  drv 


cv,  =  _ NO 

rm  fcrV 


*$572  I  7=^-  /  [(X“UT" 

1  0  v  v  0 


w _ -NO  r  r  (rv> 


Cp/  = 


■  *****  Jo  ^  ■{) 


I  j=-  /[%  +{x-{Ur  )QmN+i/2^2) 


+  §  ~  ix+{Ur^-3/2^  ~  8?  Q^-l/2^2)]e“<mNnT  dT  dr 


(12) 
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where 


a>„  =  1  + 


(x-Urr+(r-r  )c 


Srr, 


(13) 


For  r  >  R  ,  the  coefficients  for  the  radial  components  may 
he  identified  with  previous  results  [9]. 

In  general,  further  analytical  simplification  of  these 
coefficients  to  a  simple,  closed  form  involving  no  integra¬ 
tions  is  not  possible.  Still  they  represent  an  essential 
simplification.  That  is,  for  any  harmonic  m  only  Q^+^/2  * 
QmN-l/2  30(1  ®mN-3/2  are  re^uired*  or  equivalently 

W  '  Vlfl  '  °mN-3/2  “d  W  £t°"  MCMrr0Bee 

formulas  (10).  This  is  in  contrast  to  results  in  termB  of 
elliptic  integrals  [11]  in  which  the  number  of  necessary  func¬ 
tions  to  be  computed  rises  in  geometric  progression  with  m  . 

The  Legendre  function  of  second  kind  and  half  order  has 
been  examined  in  detail  [10]  in  connjction  with  recent  work  on 
the  ducted  propeller  and  extensive  tables  for  orders  of  -1/2 
through  21/2  have  been  computed  [12J.  For  a  fixed  value  of 
n  ,  has  a  s^jnPl°  logarithmic  singularity  at  a>  “  1 

and  decays  monotonically  to  zero  as  uj  increases  to  infinity. 
In  addition,  Qn  decays  monotonically  to  zero  with 

increasing  n  for  a  given  cd  • 

Altogether,  it  appears  that  the  Legendre  function  of 
second  kind  and  half-integer  order  is  the  "natural"  function 
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for  the  study  of  the  velocity  field  of  a  propeller. 

STEADY  VELOCITY  FIELD 

If  we  examine  that  part  of  the  velocity  field  correspond¬ 
ing  to  ra  =  0  ,  considerable  simplification  in  the  complex 
Fourier  velocity  coefficients  is  possible.  This  part  corre¬ 
sponds  to  the  steady,  or  time-independent  part  of  the  velocity 
field.  Using  Up  ...  to  designate  these  steady  velocity 
components,  we  obtain  from  Eqe.  7.  9,  10  and  12  together  with 
the  identity  Qjygta)  =  0^3/2 (“)  '  that 


ur  =  0 


vr  =  0 


w  =  -JSS- 

r  47T2- 


*372  f  7375- Q^H>  drv 


0  V 


H  ,  t  * 

r  (O 


Up  '  = 

r  V 


072-  f  ~7=r-  f  (rQ^m2)  "  rvQV®2)J  dT  drv 

i  A  V  A 


0  V*v  0 


R 


f  f  (x"Ut)  Qj*(a)2)  dT  drv 


Nn  J  •  '-v 

V'  ^^75  J 


0  v  v  0 


372  /  “7172  f  fVW  "  1Q-55(q)2)]  dT  drv 


Wp,  = 

r  ii_c 


V 


0  v  0 


(14) 
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As  a  result,  we  find  that  the  hound  blade  vortices  contribute 
to  the  tangential  velocity  only. 

The  integrations  over  r  occurriv.-j  in  the  expressions  for 
Up-  ,  Vp ,  and  Wp/  can  be  carried  out  analytically.  Since  it 
is  the  simplest,  we  first  consider  Vp,  and  note  the 
relationship 

(x-Ut)  Q^2)  =  -  ^  Q^(fflg)  (15) 

to  obtain 

R 

v  =  47w?  {  r'(Tv)  ^  V»i>  drv  (16) 

Since  x  appears  only  as  a  squared  quantity  in  the  argument 
cop  given  by  Eq.  11,  we  see  immediately  that  Vp,  is  symmetric 
in  x  .  From  appropriate  expansions  of  for  small  and 

large  arguments  fl2),  it  can  be  shown  that  Vp,  is  logarith¬ 
mically  infinite  over  the  propeller  disk  and  vanishes  on  the 
propeller  axis.  Because  Vp  =  0  ,  the  same  conclusions  are 
true  for  the  total  steady  induced  radial  velocity 

V  »  (Vp+Vp,)  . 

The  Up,  and  Wp,  components  are  reduced  by  rewriting 
Q^(od)  and  in  their  integral  form  [10],  or 

_ cos  2nct _  . 

[2(o)-l)+4sin2a]3/2 


°n-hM 


(17) 
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The  order  of  the  a  and  t  integrations  is  then  reversed, 
after  which  the  integrations  over  t  and  then  over  a  may  be 
carried  out.  The  final  results  are 
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Eqs.  18  can  bci  still  further  simplified  if  we  integrate 
by  parts  with  r'drv  as  the  differential  and  use  the  relation¬ 
ship  which  we  have  found, 


3r^  "  i  n  [rvQJs(®l)  +  rf,4(ml)]  ; 


+  IF  r<r , 


-  IF  r>r. , 


Since  the  circulation  must  vanish  at  'she  tip,  we  find  for  Up- 


A 

Kfl  f 

^  i 


F(0 


i  ^ 


Q-SS<“1)  ?rv  ’ 


IF  r>R  ,  -oo<X<co 


o.R  r<R  ,  x<0 


V  ■  +  ^57?  I  dI'-  ’  IF  ’  300 


t  IF  r<R  ,  x=0  (22) 

From  the  first  oZ  Eqs.  14  and  Eqs.  22,  we  see  the  following. 

At  radial  stations  greater  than  the  propeller  radius,  the  total 
steady  axial  induced  velocity  u  =>  (iip+Up')  i0  antisymmetric 
with  respect  to  x  .  At  the  propeller  plane  it  vanishes  off 
the  propeller  disk  and  has  the  same  shape  as  the  blade  circu¬ 
lation  distribution  over  the  propeller  disk,  cf.  Refs.  [14]  and 
(15)  .  Far  down  in  the  slipstream,  u  is  twice  that  at  the 
propeller  plane,  the  same  as  in  momentum  theory. 
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For  Wp-  we  get 


IF  r>R  ,  -«o<x<co 
OR  r<R  ,  x<0 


NT(r) 


J 


IF  r<R  ,  x=0 


(23) 


The  total  steady  induced  tangential  velocity  w  =>  (Wp+Wp') 
from  the  third  of  Eqs.  14  and  Eqs.  23  is,  then, 


0 


r>R  ,  -«o<x<t» 
r<R  ,  x<0 


w 


Nr(r) 

2rrr 


IF  r<R  ,  x>0 


Nf(rl 

4Trr 


IF  r<R  ,  x=0 


(24) 


That  is,  the  tangential  velocity  vanishes  everywhere  outside 
the  propeller  slipstream  and  is  proportional  to  r(r)/r  inside 
the  slipstream.  This  is  the  same  result  as  we  can  derive 
directly  by  application  of  Kelvin's  Theorem  using  a  circular 
path  about  the  x-axis . 
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The  results  of  Eqs.  22  and  24  as  regards  to  the  inflow  at 
the  propeller  plane  are  especially  noteworthy.  In  particular, 
u/w  =  fir/u  and  30  from  the  velocity  diagram  at  a  blade  section, 
the  resultant  induced  velocity  vector  is  perpendicular  to  the 
resultant  free  stream  vector  (Ui^-firi^)  .  This  is  the  same 
as  Moriya  [73  found  for  the  total  inflow,  or  steady  pluB 
higher  harmonics.  Therefore,  we  conclude  it  must  be  true  for 
these  harmonics  as  well. 

In  summary,  the  total  axial  component  of  the  steady  in¬ 
duced  velocity  field  is  determined  by  Eqs.  22  and  Up  «*  0  ,  the 
radial  component  by  Eq.  16  and  v„  =  0  ,  and  the  tangential 
component  by  Eqs.  24.  We  see  that  these  relations  if  expressed 
in  non-dimensional  form  would  be  Independent  of  the  blade 
number  for  a  fixed  advance  ratio  j  =  u/flR  and  disk  loading 


dc 

_  T  o  2N  r  _T_ 
cl  ( r/R )  7TJ  R  UR 


(25) 


whore  CT  is  the  propeller  thrust  coefficient.  Consequently, 
on  a  physical  basis  they  must  be  equal  to  the  respective  com¬ 
ponents  of  the  (Sp+Sp")  of  EQ3'  7  subject  to  the  same 

conditions.  This  limit  corresponds  to  the  velocity  field 
associated  with  what  we  now  define  as  the  GENERALIZED  ACTUATOR 
DISK,  i.  e.,  the  precise  mathematical  definition  of  the  actua¬ 
tor  disk  as  opposed  to  a  model  with  certain  characteristics 
assumed  a  priori. 
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REPRESENTATIVE  BLADE  CIRCULATION  DISTRIBUTION 
VELOCITY  PROFILES 

To  illustrate  the  theory,  we  have  calculated  several  u 
and  v  profiles  for  the  representative  blade  circulation  dis¬ 
tribution  given  by 


This  particular  r  was  contrived  because  it  simplifies  the 
calculations  and  approximates  the  familiar  Goldstein  optimum 
distribution  [6]  quite  well,  particularly  the  proper  square 
root  behavior  at  the  tip,  see  Fig.  2.  The  constant  A  is 
proportional  to  the  propeller  thrust  coefficient  cT  ..  By 
substituting  Eq.  26  into  Eq.  25  and  integrating,  we  have 


A 


7TJ 

IT 


(27) 


O  O 

where  CT  =  T/ttR  JjpU  .  T  is  the  propeller  thrust  and  p  is 
the  fluid  density. 

The  axial  and  radial  velocities  were  calculated  on  a 
CDC  1604  digital  computer  from  Eqs.  22  and  Eq.  16.  The  results 
in  terms  of  cT  from  Eq.  27  are  tabulated  in  Tables  1  and  2 
and  sketched  in  Figs.  3  and  4  respectively,  several  features 
arc  prominent  in  addition  to  those  observed  previously.  From 
Fig.  3  the  steady  induced  axial  velocity  off  the  propeller  disk 
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COMPARISON  OP  STEADY  RADIAL  VELOCITY  PROPILES 


vanishes  in  the  propeller  plane,  leaving  the  free  stream  unper¬ 
turbed  except  for  the  radial  velocity,  inside  the  slipstream 
5  quickly  approaches  within  a  blade  radius  twice  its  value  at 
the  propeller  plane  and  decays  just  as  fast  upstream.  Prom 
Fig.  4  we  see  a  rapid  change  in  the  radial  velocity  in  the 
immediate  vicinity  of  the  propeller  plane  and  a  smoother  vari¬ 
ation  farther  away.  Within  a  certain  radius  from  the  propeller 
axis  and  near  the  propeller  plane,  v  reverses  and  the  radial 
flow  is  outward.  This  reversal  is  due  to  the  influence  of  the 
trailing  vortices  inboard  of  the  maximum  value  of  r  which 
are  of  opposite  sign  to  those  outboard.  At  axial  stations  near 
the  propeller  plane  the  velocity  is  still  appreciable  for  large 
radial  distances,  whereas  the  profile  quickly  decays  with  large 
axial  distance  away  from  the  propeller  in  the  vicinity  of  the 
propeller  axis. 

CONSTANT  BLADE  CIRCULATION  DISTRIBUTION 
VELOCITY  PROFILES 

For  the  special  case  in  which  the  circulation  distribution 
over  the  propeller  blade  is  constant,  the  results  for  tho 
steady  induced  velocities  can  be  put  in  simple,  closed  form 
because  fundamentally  they  are  the  "building  blocks"  for  the 
case  of  arbitrary  circulation. 

With  T(r)  =  r  ,  it  is  easiest  to  reduce  u  by  returning 
c 

to  the  intermediate  form  of  Up/  or  the  first  of  Eqs.  18. 


This  yields  only  a  contribution  from  the  tip  which,  together 
with  Up  =  0  ,  gives 

( 

IF  r<R  ,  x<0 


Nr  n  r  _ 


OR  r<R  ,  x>0 


hTen 

4tt2u 


°-!*<a>3)  -  I  VW]  ’ 


IF  r>R  ,  x<0 


OR  r>R  ,  x>0 


(28) 


where 


03  =  i  + 

3  2rR 


02  s  sin 


-1 


k„  s 


4rR 


^(r-R)2  '  2  V*2+(r+R) 


(29) 


Across  the  cylinder  (r=R  ,  x>0)  ,  there  is  a  constant  jump  in 

velocity  of  Nrcfi/27rU  Which  can  be  identified  as  the  average 

vorticity  per  unit  length  in  the  axial  direction  resulting  from 

the  trailing  tip  vortices  or  (Nrc/27rR)/J  .  Over  the  propeller 

disk  u  becomes  Nr  0/4'/rU  ,  or  one-half  of  its  value  far 

c 

downstream,  the  same  as  in  momentum  theory. 

The  radial  component  follows  from  Eq.  16  and  Vp  »  0  in 


similar  fashion,  or 
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As  opposed  to  the  representative  case,  v  now  is  finite  over 
the  propeller  disk  except  at  the  tip.  Furthermore,  it  does  not 
reverse  sign  as  the  only  trailing  vortices  of  opposite  sign  to 
the  tip  vortices  lie  along  the  propeller  axis  and  do  not  induce 
any  radial  flow. 

These  velocities  have  been  calculated  and  are  tabulated  in 
Tables  3  and  4.  They  are  compared  back  with  the  results  for 
the  representative  circulation  distribution  in  Figs.  3  and  4. 
From  the  figures  we  see  that  for  radial  stations  less  than  the 
propeller  radius  the  results  are  not  in  close  agreement,  the 
disparity  being  most  pronounced  within  the  propeller  slipstream. 
However,  there  is  no  significant  difference  between  the  two 
results  for  r/R  >1.5  . 

While  the  condition  of  constant  circulation  along  the 
ontire  blade  is  physically  impossible  [14],  we  can  conclude 
that  the  determination  of  the  propeller  induced  velocities  out¬ 
side  the  slipstream  may  be  satisfactorily  approximated  by 
assuming  such  a  circulation  distribution.  On  the  other  hand, 
prediction  of  the  flow  field  within  the  slipstream  as  well  as 
in  the  immediate  vicinity  ahoad  of  the  propeller  plane  will  be 
in  serious  error  if  a  constant  blade  circulation  is  assumed. 

IDENTIFICATION  WITH  ACTUATOR  DISK  SOLUTION 

As  wo  said  earlier,  the  assumption  is  often  made  simply 
a  priori  that  the  blade  number  N  is  infinite  and  so,  the 
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propeller  is  replaced  by  an  "imaginary11  disk  across  which  the 
axial  velocity  is  continuous  while  the  fluid  pressure  is 
suddenly  increased  in  passing  from  one  side  to  the  other,  in 
the  absence  of  swirl  this  permits  the  use  of  appropriate  dis¬ 
tributions  of  ring  vortices  or  sources  to  represent  the 
flow  [l6].  These  representations  can  be  related  to  our  solu¬ 
tion  for  the  steady  induced  flow  field  of  the  finite-bladed 
propeller. 

First  consider  the  ring  vortex  representation  for  the 
axial  component  of  the  induced  velocity.  The  stream  function 
¥  for  a  ring  vortex  of  unit  strength  is 


V  *s  — -  J  TV  O.  /(0,  ^ 

•cen  “v  ~V  4' 


®4  "  1  + 


(x-xv)2+(r-rv)2 
2rr. . 


(31) 


where  xv  is  the  axial  location  of  the  ring  and  rv  is  the 
ring  radius.  The  corresponding  induced  axial  velocity  becomes 
after  simplification, 


1  S'? 
r  §r 


S 


2*r 


^7r 


[  rQ^(o>4 ) 


(32) 


With  t  replaced  by  x ^/\i  ,  we  see  that  u  from  Eqs .  14  is 
equal  to  the  slipstream  integration  of  Eq.  32  weighted  by  the 
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vortex  strength  -N2df(rv)dxv/27rU  .  But  this  agrees  exactly 
with  the  usual  expression  for  the  induced  axial  velocity 
associated  with  the  actuator  disk  in  terms  of  ring  vortices: 
Take  for  simplicity  the  case  of  uniform  loading  which  requires 
only  a  single  semi-infinite  solenoid  of  such  vortices  of  radius 
R  .  Prom  momentum  theory,  the  ratio  of  the  total  velocity  in 
the  ultimate  slipstream  to  u  is  V 1+CT  ,  or  linearized  for 
light  loading  (1+^0^)  ;  cf.  Table  3,  u/uc^,  =  0.25  at  x“0  . 
r/R<l  and  u/UCT  -♦  0.5  as  jc-*»  ,  r/R<l  .  Consequently, 
^UCTdxv  is  the  proper  vortex  strength.  On  the  other  hand, 
dr(rv)  =■  -rc  at  R  and  the  vortex  strength  from  the  express¬ 
ion  above  is  then  simply  NflT^dx '.^/SKO  which  properly  corre¬ 
sponds  to  the  velocity  jump  noted  for  Eqs.  28  times  dxv  . 
Substituting  CT  =  (N/mJ)  (rc/lJR)  «  NOI^/mU2  from  integration 
of  Eq.  25,  we  arrive  at  l*ucTd::v  again.  The  case  of  arbitrary 
loading  follows  from  the  superposition  of  such  solenoids  of 
radius  r 

v 

If  we  examine  in  turn  the  radial  velocity  of  a  vortex  ring 
or  d'F/rdx  from  Eq.  31  and  compare  it  similarly  with  v  from 
Eqs.  lJ»,  we  obtain  the  identical  result.  With  regard  to  w  , 
the  swirl  is  generally  omitted  from  the  actuator  disk  concept. 
But  if  we  incorporate  a  distribution  of  concentric  cylinders  of 
semi-infinite  straight  vortices  parallel  to  the  x-axis,  we  can 
establish  the  equivalence  as  for  u  and  v  . 
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For  the  representation  of  the  actuator  disk  by  the  distri¬ 
bution  of  ring  sources  on  the  propeller  disk,  we  can  proceed 
formally  to  prove  the  equivalence  as  for  the  vortex  rir.g3.  In¬ 
asmuch  as  the  equivalence  between  the  ring  sources  and  vortex 
rings  has  already  been  established  [ 16]  &  [17],  this  is  not 
necessary.  If  desired,  though,  Eqs.  22  are  best  for  the 
axial  velocity.  For  the  radial  velocity,  integrate  Eq.  16  by 
parts.  The  respective  velocities  for  a  ring  source  follow  from 
the  potential  <I>  per  unit  source  strength, 

<t>  e  _  -i—  ■  *  q  ,  (a>. )  (33) 

4 IT  </rr^  1 

where  x  =  0  is  the  axial  location  of  the  ring  and  rv  ,  the 
radius,  we  find  that  outBide  the  propeller  slipstream  the 
velocities  are  equal  and  the  required  strength  of  the  ring 
source  is  -Nnrvr(rv)dry/U  ,  i.  e.,  a  ring  sink.  Inside  the 
slipstream  the  radial  velocities  are  the  same  but  the  axial 
velocities  differ  by  the  constant  NfiT(r)/2?iU  ,  cf.  Eq.  22, 
since  u  for  the  ring  source  is  antisymmetric  in  x  every¬ 
where  . 


CONCLUSIONS 

From  our  study  of  the  induced  velocity  field  of  a  finlto- 
bladed  propeller  with  arbitrary  circulation  distribution,  we 
have  concluded  the  followings 
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Relatively  simple  forms  for  the  Fourier  coefficients  of 
the  axial,  radial  and  tangential  velocities  can  be  obtained 
involving  only  Legendre  functions  of  the  second  kind  and  half 
integer  order.  These  appear  to  be  the  "natural"  functions  for 
use  in  propeller  theory. 

A  comparison  of  the  steady  velocity  profiles  for  both  a 
representative  and  a  constant  circulation  distribution  reveals 
that  the  two  results  in  general  agree  closely  outside  the  pro¬ 
peller  slipstream.  However,  calculation  of  the  field  immedi¬ 
ately  ahead  of  the  propeller  as  well  as  inside  the  slipstream 
under  the  assumption  of  uniform  loading  may  lead  to  serious 
errors . 

For  the  special  case  of  constant  circulation,  the  steady 
induced  velocity  components  cm  be  expressed  in  closed  form 
involving  no  integration. 

The  axial  and  radial  induced  velocities  for  the  conven¬ 
tional  actuator  disk  are  established  as  exactly  equal  to  their 
steady  counterparts  for  a  finite-bladed  propeller,  provided 
the  advance  ratio  and  the  disk  loading  are.  the  same. 
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PRINCIPAL  NOMENCLATURE 

p  p 

propeller  thrust  coefficient,  T/rrR  ?<pu 

complex  Fourier  velocity  coefficients  for 
gp  and  2p; 

propeller  advance  ratio,  u/ftR 

harmonic  number 

number  of  propeller  blades 

Legendre  t unction  of  second  kind  and  half- 
integer  order  of  argument  cd 

vector  velocity  induced  by  bound  blade  vortices 

vector  velocity  induced  by  blade  trailing 
vortices 

propeller  radius 

forward  flight  velocity 

axial,  radial  and  tangential  components 
respectively  of  induced  velocity 

cylindrical  propeller-fixed  coordinate  system 

propeller  blade  circulation 

Heuman's  Lambda  function  of  argument  0 
and  modulus  k 


propeller  angular  velocity 


differentiation  of  a  function  with  respect 
to  its  indicated  argument 


zeroth  harmonic  or  steady  part  of  induced 
velocity 


Approved  Distribution  List  for  Unclassified 
Technical  Reports  Issued  Under  Contract  Nonr-4357 (00) 


Chief,  Bureau  of  Naval  Weapons 
(RAAD-3) 

Department  of  the  Navy 
Washington,  D.  C.  20360 

Chief,  Bureau  of  Naval  Weapons 
(RAAD-22) 

Department  of  the  Navy 
Washington,  D.  C.  20360 

Chief,  Bureau  of  Naval  Weapons 
(RAAD-32) 

Department  of  the  Navy 
Washington,  D.  C.  20360 

Chief,  Bureau  of  Naval  Weapons 
(RAAD-33) 

Department  of  the  Navy 
Washington,  D.  C.  20360 

Chief,  Bureau  of  Naval  Weapons 

(««««««  **  t  \ 

Dona  rf  a  f 

- wmwttM  W4.  MMTJ 

Washington,  D.  C.  20360 

Chief,  Bureau  of  Naval  Weapons 
(RA-4) 

Department  of  the  Navy 
Washington,  D.  C.  20360 

Chief,  Bureau  of  Naval  Weapons 
(R-55) 

Department  of  the  Navy 
Washington,  D,  C.  20360 

Chief,  Bureau  of  Naval  Weapons 
•KRRE-4) 

Department  of  the  Navy 
Washington,  D.  C.  20360 

Commanding  Officer 
U.  S.  Army  Transportation  Research 
Command 

Fortr  East  is,  Virginia 
ATTN:  SMOFE-TD  (1  copy) 

Research  Reference  Center  (1 

Chief  of  Naval  Research  (Code  46.1) 
Department  of  the  Navy 
Washington,  D.  C.  20360  (6  copies) 

Chief  of  Naval  Research  (Code  438) 
Department  of  the  Navy 
Washington,  D.  C.  20360 


Commanding  Officer 
Office  of  Naval  Research  Branch 
Office 

Navy  #100,  Eox  39,  F.P.O. 

New  York,  New  York  (2  copies) 

Commanding  Officer 
Office  of  Naval  Research  Branch 
Office 
346  Broadway 
New  York  13,  New  York 

Director 

Naval  Research  Laboratory 
Technical  Information  Office 
Washington,  D.  C.  20390  (6  copies) 

Commander 

Army  Material  Command 
ATTN:  AMCRD-RS-PE-A 
Department  of  the  Army 
Washington,  D.  C.  20315 

Commanding  Officer  and  Director 
David  Taylor  Model  Basin 
Aerodynamics  Laboratory  Library 
Washington,  D.  C.  20007 

Army  Research  Office 
Physical  Sciences  Division 
ATTN:  Mr.  R.  Ballard 
3045  Columbia  Pike 
Arlington,  Virginia  20310 

U.  S.  Air  Force  (SRGL) 

Office  of  Scientific  Research 
Washington  25,  D.  C. 

Research  and  Technology  Division 
ATTN:  SESSC  (Mr.  Lindenbaum) 
Wright-Patterson  AFB,  Ohio  45433 

Defense  Documentation  Center 

Us*  <1 

•*^»  /  WWtllUkVIl  M  Vt* 

Building  #5 

copy)  Alexandria,  Virginia  22314 
(20  copies) 

National  Aeronautics  and  Space 
Administration 

600  Independence  Avenue,  S.W. 

ATTN:  Code  RA,  Code  RAD 
Washington,  D.  C.  20546  (2  copies) 


Library 

American  Institute  of  Aeronautics 
and  Astronautics 
Two  East  64  Street 
New  York  21,  New  York  (2  copies) 

Research  and  Technology  Division 
ATTN:  EDM  (Mr.  Antonatos) 
Wright-Patterson  APB,  Ohio  45433 

Cornell  Aeronautical  Lab.,  Inc. 

4455  Genesee  Street 
Buffalo  21,  New  York 
atypmj  Technical  Director 
Mr.  H.  A.  Chellek 

Collins  Radio  Company 
Cedar  Rapids,  Iowa 
ATTN:  Dr.  A.  Lippisch 

Mr.  A.  M.  o.  Smith,  Supervisor 
Aerodynamics  Research  Group 
Aircraft  Division 
Douglas  Aircraft  Co.,  Inc. 

Long  Beach,  California 

Georgia  institute  of  Technology 
School  of  Aerospace  Engineering 
Atlanta  13,  Georgia 
ATTN:  Mr.  D.  W.  Dutton 

Vehicle  Research  Corporation 
1661  Lombardy  Road 
Pasadena,  California 
ATTN:  Dr.  Scott  Rethorst 

University  of  Virginia 
Aerospace  Engineering  Department 
Charlottnville,  Virginia 
ATTN:  Dr.  G.  B.  Matthews 

Mississippi  State  University 
Engineering  and  Industrial  Research 
Station 

State  College,  Mississippi 
ATTN:  Dr.  J.  J.  Cornish 

Vidya  Division 
1450  Page  Mill  Road 
Stanford  Industrial  Park 
Palo  Alto,  California 
ATTN:  Dr.  J.  N.  Nielsen 


Syracuse  University 
Mechanical  Engineering  Department 
Syracuse,  New  York 
ATTN:  Dr.  S.  Eskinazi 

Naval  Postgraduate  School 
Aeronautical  Engineering 
Department 

Monterey,  California 
ATTN:  Dr.  R,  Head 

Technical  Library  (Code  P80962) 

U.  S.  Naval  Ordnance  Test  Station 

Pasadena  Annex 

3202  E.  Foothill  Blvd. 

Pasadena  8,  California 

Commanding  Officer  and  Director 
David  Taylor  Model  Basin 
Aerodynamics  Laboratory 
Washington,  D.  C.  20007 
ATTN:  Mr.  H.  Chaplin 

Mr.  Maurice  Sevik 
Ordnance  Research  Laboratory 
Pennsylvania  State  University 
P.  0.  Box  30 

University  Park,  Pennsylvania 

Mr.  Paul  Granville 
Hydrodynamics  Laboratory 
David  Taylor  Model  Basin 
Code  581 

Washington,  D.  C.  20007 

Commanding  Officer  and  Director 
David  Taylor  Model  Basin 
Hydrodynamics  Laboratory 
Washington,  D.  C.  20007 
ATTN:  Dr.  William  B.  Morgan 


2 


